Resonance-like nuclear processes in solids: 3rd and 4th order processes 
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It is recognized that in the family of heavy charged particle and electron assisted double nuclear 
processes resonance-like behavior can appear. The transition rates of the heavy particle assisted 3rd- 
order and electron assisted 4th-order resonance like double nuclear processes are determined. The 
power of low energy nuclear reactions in Ni — H systems formed in Ni placed in H2 gas environment 
is treated. Nuclear power produced by quasi-resonant electron assisted double nuclear processes in 
these Ni — H systems is calculated. The power obtained tallies with experiments and its magnitude 
is considerable for practical applications. 
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I. INTRODUCTION 

Since the "cold fusion" publication by Fleischmann 
and Pons in 1989 a new field of experimental physics 
has emerged. Although even the possibility of the phe- 
nomenon of nuclear fusion at low energies is doubted by 
many representatives of mainstream physics, the quest 
for low-energy nuclear reactions (LENR) flourished and 
hundreds of publications (mostly experimental) have 
been devoted to various aspects of the problem. (For 
the summary of experimental observations, the theoreti- 
cal efforts, and background events see e.g. 0,0].) The 
main reasons for revulsion against the topic have been: 
(a) according to standard knowledge of nuclear physics 
due to the Coulomb repulsion no nuclear reaction should 
take place at energies corresponding to room tempera- 
ture, (b) the observed extra heat attributed to nuclear 
reactions is not accompanied by the nuclear end prod- 
ucts expected from hot fusion experiences, (c) nuclear 
transmutations were also observed, that considering the 
Coulomb interaction is an even more inexplicable fact at 
these energies. 

The situation is further complicated by the fact that 
the electrolysis, gas discharge and/or high pressure gas 
environments that are stipulated to induce LENR have 
their effect through rather complex microscopic processes 
that in most of the cases are difficult to reproduce. As 
a result, for explaining the riddle of cold fusion it is in- 
dispensable to understand theoretically the underlying 
nuclear reactions. 

In our opinion we made progress in the theoret- 
ical clarification of the nuclear physics behind LENR. 
The idea is based on the fact that if a heavy, charged 
particle (proton, deuteron) of low energy enters a solid 
(metal) , then through the Coulomb interaction it changes 



the state of charged particles, primarily quasi free valence 
electrons in the metal while its own state is also changed. 
Our results of standard perturbation theory calculations 
indicate that by means of the Coulomb interaction the 
ingoing charged particle in this change of state can ob- 
tain so high value of virtual momentum (energy) that is 
enough to induce various nuclear reactions including fu- 
sions and/or transmutations. With the results of our the- 
oretical work, we found reasonable answers to the above 
questions (a, b and c). 

In however, it was not taken into account that the 
emergence of new charged particles may alter the state 
of the solid too in a way that can allow further nuclear 
processes, such as the family of double nuclear processes 
which has special interest from the point of view of low 
energy nuclear processes. In this paper this effect is con- 
sidered. First a 3rd-order double nuclear process is inves- 
tigated the graphs of which can be seen in FIG. 1. It is 
shown that a resonance can appear in this process. Next 
a 4th-order double nuclear process is discussed the graphs 
of which can be seen in FIG. 2 which shows resonance- 
like characters too. Finally the nuclear power produced 
in Ni — H systems is calculated. 

This paper is mainly based on 0. The notation 
and the outline of the calculation is the following. The 
Coulomb coupling strength e 2 = afhc and the strong 
coupling strength f 2 — 0.08/ic 0, e is the elementary 
charge, a/ is the fine structure constant, h is the reduced 
Planck constant and c is the velocity of light. 

When calculating the matrix elements of the strong in- 
teraction potential, the long wavelength approximation 
|99(0)| = fjk(k)/VV of the Coulomb solution ip(r) is 
used, that is valid in the range of a nucleon, where fjk (k) 
is the appropriate Coulomb factor corresponding to the 
particles, which take part in strong interaction and V is 
the volume of normalization. We introduced the follow- 
ing notation 
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F jk (k) - f 2 k (k) - exp[27TVjk{k)] _ 1 
Here the Sommerfeld parameter rjjk (k) for particles j 
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and k of electric charge numbers zj and Zk is determined 

as 

rj jk (k) = ZjZkaf-^-, (2) 

where fc = |kj- — kfc | is the magnitude of the relative wave 
vector k = kj — k^ of the interacting particles of wave 
vectors kj and k^ and Hjf. = mjin^/ (rrij + m&) is the 
reduced mass of particles of rest masses rrij and m&. 

For quasi-free particles (electron and ingoing proton) 
taking part in Coulomb interaction we use plane waves. 
Thus their Coulomb matrix element is calculated in the 
Born approximation which is corrected with the so called 
Sommerfeld factor 

9S (kinj k ou t) = fjk(kin) I fjk(kout) (3) 

@ , where ki n and fc oll * are the magnitudes of the relative 
wave numbers before and after Coulomb scattering. For 
other details and notation see Pf. 



II. RESONANCE-LIKE HEAVY CHARGED 
PARTICLE AND ELECTRON ASSISTED 
DOUBLE NUCLEAR PROCESSES 

Preliminarily one must emphasize that in both (3rd 
and 4th order) double nuclear processes discussed res- 
onances may occur. The reason for the possibility of 
resonance is that the continuum of the kinetic energy of 
an intermediate state is shifted down by the energy of 
the (first) nuclear transition and therefore one of the de- 
nominators in the perturbation calculation can be equal 
to zero. The occurrence of resonances increases the rate 
significantly. 




FIG. 1: The graphs of quasi-resonant heavy charged parti- 
cle Coulomb-assisted nuclear reactions. The filled dot de- 
notes Coulomb-interaction and the open circle denotes nuclear 
(strong) interaction. Free particle 2 (e.g. p or d) belongs to 
system B. Localized particles 1, 3 and 4 belong to system A. 
All the particles are heavy and positively charged. From the 
point of view of the nuclear process particles 1' and 2' are 
ingoing, particles 3 and 4 are targets and particles 5 and 6 
are reaction products. In the case discussed in the text in the 
process of FIG. 1(a) there is a resonance, therefore its con- 
tribution to the rate is dominant and the contribution of the 
process of FIG. 1(b) is neglected. 



A. Resonance-like heavy charged particle assisted 
double nuclear processes 

First the resonance- like heavy charged particle assisted 
double nuclear processes (see FIG. 1) are discussed. The 
particles in FIG. 1 are all heavy, and positively charged. 
The ingoing particle is particle 2, which belongs to system 
B (the ensemble of incoming particles forms system B). 
It is supposed that it has moderately low energy (of about 
keV order of magnitude), that is raised in the second or- 
der processes discussed in [ij or after it in the deceler- 
ating process. Particle 2 scatters by Coulomb scattering 
on particle 1 localized in the solid (system A). Particles 
3 and 4 are the nuclear targets of system A. The parti- 
cles in the intermediate state 1' and 2' may pick up large 
enough wave vector to overcome the Coulomb repulsion 
due to particles 4 and 3. Particles 5 and 6 are products of 
the process. The nuclear process 1' + 4 — > 5 takes place 
as the consequence of the modification of system A by 
system B. Here both nuclear processes are thought to 
be nuclear captures. It can be shown (see Appendix I.) 
that the process may have resonance like character if the 



masses of particles 5 and 6 differ significantly, therefore 
the contribution of the leading graph (e.g. FIG. 1(a) in 
the case m 5 ^> m 6 discussed) to the rate is enough to 
calculate. 

We take as initial state of a localized particle k 

i> kl (x) = (Pl/nf 4 exp (-/?, 2 x 2 /2) , (4) 

which is the ground state of a 3-dimensional harmonic 
oscillator of energy — ^fkuk and of angular frequency 
Wfe 0- Now we take particle 1 (k = 1) localized and the 
initial state ipu (x) is of the form ([3]) with the parameter 
Pi = y/m\iij\/h, where mi is the rest mass of the lo- 
calized particle 1. The Coulomb matrix elements which 
contain initial state of form ip^i (x) also preserve momen- 
tum Q. 

The total rate of the reaction 

W^l = K^K[ 3) (A,A 5 )h 2 corr ^h 2 corrA N 13i (5) 
xJ2H(k 2t ,k r )N 2 (k 2l ), 

k 2 i 
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where N134 is the number of triples of particles 1, 3 and 
4, that may be effective for one particle 2, and N 2 (k 2 i) 
is the actual number of quasi-free particles 2 of wave 
number k 2 i- Furthermore 



= 92l6a 2 f n 7 / 2 [ 1 



He 



R 2 cpl (6) 



where R = 1.2 x 10 [cm] is the proton-radius, 



K f) (A> As) = ^kSWA_^6 ; (?) 



V2 A 5 -^A 



m 2 



where A5 = A01 + A04 — A05 is the energy of reaction 
1' + 4 — > 5, and A 6 = A 02 + A 03 — A 06 is the energy of 
reaction 2' + 3 — > 6. The Aoj-s are the energy defects of 
the corresponding nuclei and A = A 5 + A 6 is the total 
reaction energy. 



H{k 2h k r )=g 2 s {k 2U 2k r )F^{k r )F li {k r l (8) 



where 



k r 



2m 2 (A 5 -^A 



h 



(9) 



/l, orr 3, tlrnrr 4 are defined by Eq.(45) of Q and are de- 
termined as 



hcorr,k A]^ 

in the case of proton captures, and 



^corr,k 



Ai 



(10) 



(11) 



in the case of dcutcron captures. (For the details of the 
calculation see Appendix. I.) 

In a numerical example let particles 1, 2, 3 be 
deuterons, particle 4 be Iq 6 PcI, particle 5 be l^ 8 Ag and 
particle 6 be 2 He. If particle 1 is a deuteron then 
/3i = 4.81 x 10 8 [cm- 1 ] (see Sec. VIII. A of [4]). With 
these choices 

A 6 = 23.847 [MeV] , A 5 = 10.833 [MeV] , (12) 



A = 34.681 [MeV] 



(13) 



^ 12 c 2 = m Q c 2 = 931.494 [MeV] , /i 56 c 2 ~ 4m c 2 , (14) 



Vm/ma = 1/28, ^ 56 /m 2 ~ 2. 
With these numbers one can obtain 

^56 



A 5 



m 5 



-A = 9.60 [MeV] , 



K[ 3) (A, A 5 ) = 5.00 x 10 5 



(15) 

(16) 
(17) 



and 

hck r = 189 [MeV] , (18) 
the latter producing 

F 23 (k r ) = 0.89, F 14 {k r ) = 2.79 x 10~ 8 . (19) 
Finally one gets 

W^ t = 1.45 x 10 8 x iV 134 9s ( k ^ 2kr ) N 2 (*2i) t^ 1 ] • 

(20) 

The (k 2 i,2k r ) dependence is similar to the 
Gs(k 2 i,2kif) dependence discussed at the end of 
section VIII. B of [J. For the number iVi34, in the cases 
discussed, as a lowest estimation 



^134 = (V e ff/v c )r Al r A3 r A4 , 



(21) 



where r Ai is the relative natural abundance of isotopes 
of mass number Ai or r Ai = u if the z-th particle is a 
hydrogen isotope, V e // is the volume effectively felt by 
a particle 2, u is the proton (or deuteron) over metal 
number density and v c (= d 3 /4) is the unit cell of the 
solid (in the case of fee metals). 

It should be noted that there may be a great variety 
of other types of possible quasi-resonant heavy charged 
particle assisted double nuclear reactions the discussion 
of which is not given here. 



B. Resonance-like electron assisted double nuclear 
processes 

From the processes discussed up till now one can con- 
clude the following: (a) it is advantageous, if the wave 
vector (momentum) transferred through the intermediate 
state by Coulomb interaction has the possible maximum 
value, (b) the electron assisted process is advantageous 
since the Coulomb and Sommerfeld factors cause minimal 
or negligible hindering in this case and (c) the appear- 
ance of resonance significantly increases the rate. These 
conclusions led us to a 4th-order, electron assisted dou- 
bled nuclear process, the graphs of which can be seen in 
FIG. 2. At a particular choice of the participants (see 
Appendix II.) resonances can be found in the processes 
of FIG. 2 (a) and (b). The resonance arises in line 4'. 
Moreover, the details of the calculation show that high 
contribution to the rate is obtained if the energy of the 
electron is negligible in the final state compared with the 
energy of the total reaction. In other words, the main 
contribution to the rate is produced by final states in 
which particles 6 and 7 share the reaction energy. 

It is supposed that particle 1 (1' and 1") is a quasi-free 
electron of the solid (a metal). Particle 2 is an ingoing 
free, heavy, positively charged particle of system B. Par- 
ticles 3, 4, and 5 are heavy particles of positive charge 
that are localized in the solid. If processes of FIG. 2 (a) 
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FIG. 2: The graphs of quasi-resonant, electron assisted nu- 
clear processes. The simple lines represent free (initial (1), 
intermediate (1') and final (1")) electrons. The doubled lines 
represent heavy, charged initial (free 2, bound 4) particles 
(such as p, d), their intermediate states (2', 4'), target nu- 
clei (3,5) and reaction products (6,7). The filled dot de- 
notes Coulomb- interaction and the open circle denotes nuclear 
(strong) interaction. 



5, 

Kjf> = 2-3^/^) (£) 4 (l - I)" (^- 2 - 4 ) 2 (23) 
x (he) R 2 ckf Max f3l, 

where ki,Max is the maximum of the possible wave vec- 
tors of the electron in the final state. (We calculate the 
rate of those processes in which the kinetic energy E\ / 
of the electron can be neglected in the energy Ef of the 
final state). 



(Fg 



F 12 (2k r 4 ) 



F 23 (k r 4 )F 45 (fcj) 



(24) 



with 



v /2TO 4 <5 a ,h 



(25) 



Now Ag = A02 + A03 — Ao6 is the energy of reaction 
2' + 3 — s- 6 and A7 = A 5 + A 4 — A 7 is the energy 
of reaction 4' + 5 — >• 7. The Aqj-s are again the energy 
defects of the corresponding nuclei and the total reaction 
energy 



A 7 . 



X(A) 



mj v / 2m7C 2 A / 1 

<5?,A?„ I a: 



1 

A~b 



(26) 



(27) 



where 8 a ,b is given by (jlOip . furthermore, for A Q , Aj, and 
A a b and see (|104l) . (|105[) and (|106l) . and for h CO rr,3 and 
h CO rr,5 see (p~0|) and (ITT|) . 
It is reasonable to take 

^345 1 4 

—y- = —rA 3 rA 4 rA 5 = -^rA s rA A rA 5 , (28) 

where v c — d 3 /4 is the unit cell of the solid in the case 
of fee metals, r Ai is the relative natural abundance of 
isotopes of mass number Aj or r Ai = u if the i-th particle 
is a hydrogen isotope. Thus 



and (b) have resonance-like character (see Appendix II.) 
then the process of FIG. 2 (c) has not and therefore its 
contribution may be neglected. Now particle 4 is local- 
ized and its state is ipa (x) given by (j3|. 



The total rate W$ 
can be obtained from the rate wif'isee (|107[) ) as 



(for the details see Appendix II.) 

fi 



W, 



(4) 



(F G ) X (A)h 2 h 2 N 1 N, 



N345 
V : 



(22) 



where iVi is the instantaneous number of the quasi-free 
electrons that are felt by one ingoing particle 2, N 2 is the 
number of quais-free ingoing particles 2, -^345/!^ is the 
number density of the target triples of particles 3, 4 and 



Wtot = (^<T) (Fg) X (A) h I corr ^h z corr ^N 1 N 2 r A3 r Al r A , 

(29) 

If particle 4 is a deuteron then (3 4 = 4.81 x 10 8 [cm 1 ] . 
Taking E\_Max — 0.1 [MeV], that is the maximum of the 
possible energies of the electron in the final state, one 

gets hi Max 



= 3.50 x 1(T 3 (z 2 z 2 z 4 y [MeV 



d 3 



(30) 



in the case of Pd, i.e. with d = 3.89 x 1CP 8 [cm]. 

As a numerical example let particles 2, 4 and 5 be 
deuterons, particle 3 a Pd isotope (of ^3 = 106 with 
r A 3 — 0.2733) and particle 7 2 He. With this choice 



10.833 [MeV] , A 7 = 23.847 [MeV] , 



(31) 
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A = 34.681 [MeV] , (32) 

furthermore 

6 a , b = 9.60 [MeV] , A ab = 198.7 [MeV] , (33) 

A Q = 19.2 [MeV] , A b = 179.5 [MeV] , (34) 
and mi jm^ = 2 resulting 

X (A) = 1.69 x 1CT 6 [MeV' 6 ] . (35) 



F 23 (fcj) = 2 - 79 x 10- 8 ,F 45 (kl) = 0.89,F 12 (2*J) = 1 

(36) 

and 



(Fi 2 (Ai(£i))) = 23.18 x fa 1/2 [eF] 
producing 



(F G ) = 5.76 x 10~ 7 x (E, 



-1/2 



(37) 



(38) 



where Ei is the energy of the initial conduction electron. 

— 1/2 

Averaging E i by means of the Fermi-Dirac distribu- 
tion in the Sommerfeld free electron model at T = 
yields 



Ei 1 ' 2 ) =2(E F 1/2 [eV[ 



(39) 



where Ep denotes the Fermi energy [8j. With these num- 
bers 

W$ = 8.36 x 10- 11 [e f 1/2 [eV] \ x N^u 2 [s^ 1 ] . 

(40) 

This rate produces a total power = CpW^A, 

where C E = 1-602 x 10~ 13 [J /MeV] is the energy unit 
conversion factor. With a deuteron concentration inde- 
pendent E F = 17 [eV], 

P^l = 1.13 x 10" 22 x N^u 2 [W] . (41) 

In all the charged particle assisted processes discussed 
the quasi- resonant, 4th-order electron assisted double nu- 
clear process seems to be the leading one. 



III. NUCLEAR POWER IN NI-H SYSTEMS 
DUE TO DOUBLE PROTON CAPTURE 

In this section we deal with a family of 4th-order reso- 
nant, electron-assisted double nuclear processes shown in 
FIG. 2. Our aim is to show that in Ni—H systems formed 
in hydrogen gas one of the family of processes in FIG. 2. 
may have high rate and the power generated by these nu- 
clear processes is also considerable from practical point 
of view that can be calculated in our theory. In Ni — H 
systems formed in hydrogen gas extra heat production 



was observed @, [l(| whose nuclear origin was proven 
with neutron detection [TTJ] . Since in Ni — H systems 
formed in H gas there is no light particle for significant 
nuclear effect save the natural deuteron content of hy- 
drogen, according to our theory the primary process to 
generate considerable energy is the quasi-resonant elec- 
tron assisted (double) proton capture of the Ni isotopes. 
The two proton captures 

(42) 



28 Ni +P 



which arc coupled due to the quasi-resonant electron as- 
sisted process, are investigated (see FIG. 2). Particle 1 
(1' and 1") is a quasi-free electron of the metal, particle 2 
is a quasi-free ingoing proton and particle 4 is a localized 
proton. Particles 3 and 5 are different Ni isotopes and 
they have mass numbers A3 and A5, respectively. Parti- 
cles 6 and 7 are Cu isotopes of mass numbers A 3 + 1 and 
A5 + 1, respectively. Both nuclear transitions 2', 3 — » 6 
and 4', 5 — > 7 are reactions of type (|4"2")l . The process is 
called quasi-resonant electron assisted double proton cap- 
ture process. Most of the daughter nuclei ^ 9 +1 Cu decay 
by the 



A+l 
29 



Ci 



A+l 
28 



Ni + Q 



EC 



(43) 



electron capture reaction. TABLE I. of [4[ contains the 
relevant data for reactions (1421) and (|43|. As it is dis- 
cussed above at a particular choice of the participants 
resonances can be found in the processes of FIG. 2. The 
resonance appears in line 4'. 

Now particle 4 is a localized proton. We take ipu (x) 
(see (0])) as initial state of particle 4. The parameter (3^ in 
the case of a localized proton is /3J = yjm p uji/h, where 
m p is the proton rest mass and W4 is the angular fre- 
quency of the ground state of a 3-dimensional harmonic 
oscillator of energy Eq = \huj±. In NiHnr5 the energy 
of an optical phonon huj^ = 88 [meV] [12|], that results 
0*1 = 6.51 x 10 8 [cm -1 ] which is used in the calculation. 
The Coulomb matrix elements which contain initial state 
of form tp4i (x) also preserves momentum Q ■ The num- 
ber density N345/V = Ad~ 3 rA 3 rA s u (see l[2"8j0 . where u 
denotes the proton over metal number density and d is 
the length of the Ni elementary cell (d = 3.52 x 10~ 8 
[cm] ) . 

Now we reformulate the results of section II. B. Aa 3 = 
A02 + A03 — A06 is the energy of reaction 2' + 3 — > 6 
and Aa 5 = A04 + A05 — A07 is the energy of reaction 



4' 



7. The Aqj-s are the energy defects of the 



corresponding nuclei and the total reaction energy 



Aa 3 ,a 5 = Aa 3 +A A5 . 



(44) 



It was found above that 



A (A 3 + 1) A A3 - (A + 1) A A5 

Sa ^ = A 3 + A 5 + 2 (45) 

is a crucial quantity (it was defined by (jl01[) 'l. since, if 



> then resonance appears in line 4' at 



A 3 ,A 5 



(46) 
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Substituting 

ftLr,3>&rr,B = (A* " 28 ) 2 (A* ~ 28 )" ( 47 ) 

into (122p the total rate of the 4th-order, resonance-like 
electron assisted double proton capture processes has the 
form 

W$(A 3 ,A 5 ) = ±K& (A 3 - 28) 2 (A s - 28) 2 (48) 
x (Fg) xNiN 2 r Aa r A r u. 

Here N\ is the instantaneous number of the quasi-free 
electrons that are felt by an ingoing particle 2, N 2 is the 
number of quasi-free ingoing protons that can interact 
with particle triplets 3, 4 and 5, and r As and r A& are 
the relative natural abundances of Ni isotopes of mass 
numbers A3 and A5, respectively (see TABLE I. of [4]). 

K { 4) given by flH|), where R = 1.2 x 1CT 13 [cm] is the 
proton-radius. Taking E\ t Max = 0.1 [AfeV] again one 
gets 

id~ 3 K^ = 0.0118 [MeV 6 s~ 1 ] (49) 

in the case of Ni. Since /X23 — ^45 and 22^3 = Z4Z5 = 
28 therefore / ^ (A-', : , . ~ F 45 (^.aJ in (F G ). 
Furthermore, since particle 1 is an electron therefore 
F i2{2kl AlA2 ) = 1 and (El reads 

(F G ;A 3 ,a 5 ) = (FMh)) Fi 3 (kl A ^ . (50) 

Here (Fi 2 (ki(Ei))) = 46.36 x (e~ 1/2 [eV]), where E F 

denotes the Fermi energy. If 2m p c 2 S> S, true in our 
case, then (l27l) can be written as 

(A 5 + l) 3 /y2m c2A A3 ,^ 

Here moc 2 = 931.494 [MeV] is the atomic mass unit. 

One can see from TABLE I. of Q that increases 
monotonically with the increase of A. Thus 5a 3 .a 5 > if 
A3 > A5, and the total power can be written as 

P tot = N l N 2 uP Q J2 i>A 3 ,A 5 , (52) 

A 3 >A 5 

where P = 4d~ 3 K^ ] 'C^Aq 5 '. C E = 1-602 x lO" 13 
[J/MeV] is the energy unit conversion factor, Ao = 1 
[MeV] as an order of magnitude of a typical nuclear re- 
action energy value. The use of Ao makes 

^ 3 , As = (A 3 -28) 2 (A 5 - 28) 2 r A3 r As x (53) 
xAa 3 ,a 5 (F G; a 3 ,a 5 ) Xa^a^o 

dimensionless. In our case 

P = 1.88 x 10~ 15 [W] (54) 



and 

i>A 3 ,As = 6 - 48 x 10-6 x ( e f 1/2 i eV i) ■ ( 55 ) 

A3>A 5 

In J2 As>Al> i>A 3 ,A s the significant contributions are given 

by ^61,58, "062,58, V>64,58 and 064,60, with 064,58 = 6.36 X 

10~ 6 x (^Ep 1 ^ 2 [eV]j as the leading term responsible for 

98% of the effect. A hydrogen concentration independent 
Ep = 17 [eV] is used producing 

Ptot = 2.95 x 10~ 21 x iViiVsu [W] . (56) 

Now we proceed to the determination of Ni, N 2 and u. 
N\ stands for the number of valence electrons which can 
interact with a proton (particle 2). Considering that Ni 
consists of micro crystals of linear dimension of about 
(D [fim]) 1—10 /im, it is reasonable to assume that a 
proton penetrating the material "feels" all the valence 
electrons. Thus 

Nx = ^g e , (57) 

V c 

where V g = D 3 x 10~ 12 [cm 3 ] is the volume of micro 
crystals, v c = d 3 /4 is the volume of the elementary cell 
and g e = 10 is the number of valence electrons in an 
elementary cell of Ni. From it N± — D 3 [pirn] x 9.17 x 
10 11 . N 2 denotes the number of protons which can be 
considered quasi free in respect of the process. It can be 
calculated if the number density n p in gas is multiplied 
with the metal volume V 2 , where they are taken for free. 
The volume is the product of the surface F of the sample 
and the length d of the elementary cell. The result is 

^2 = 2Fn L ^d, (58) 
PoT 

where rip — 2.69 x 10 19 [cto~ 3 ] is the number density 
in the gas of pressure po = 1 [atm] and temperature 
T = 273 [°K] . The actual pressure and temperature 
are denoted by p and T. Factor 2 follows from the fact 
that the hydrogen molecule contains 2 atoms. Here the 
catalytic process producing atomic hydrogen is not con- 
sidered, it supposed, which is a rough over estimation, 
that at the Ni surface the whole gas is atomistic. 

In Q excess heat power 44 [W] is obtained from a Ni 
rod of diameter do = 0.5 [cm] and length ho = 9 [cm] at 
p = 0.5 [atm] and T = 753 \°K~\ . At this temperature 
and pressure u ~ 7.9 x 10 -5 [13|. Utilizing this and the 
values Ni and N 2 obtained above, P tot — D 3 x 0.074 
[W] from ([56)) . In case of D ~ 8.41 [/im] the result is 
44 [W] , which considering that in this simple model a 
great number of solid state processes were neglected and 
in the nuclear processes it was only the Weisskopf ap- 
proximation in which the matrix element was calculated, 
is a very good approximation. 

The surface area of the sample, as can be seen from 
(1551) . can significantly increase the rate of the processes 
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in the sample and the power generated. For example, if 
we take a sample of similar volume to the previous one, 
but not in the shape of a rod but powder of particles of 
diameter D p then the power is increased by a factor of 
3d /(2D p ). If D p = [cm] then this factor is 750, 
that is a sample of similar volume can produce 33 [kW] 
instead of 44 [W] . It seems that this is what [l4| found 
experimentally who applies Ni powder in his patented 
device. All this means that the above mechanism can 
produce power generation for practical applications. 



IV. SUMMARY 

Resonance-like heavy particle and electron assisted 
double nuclear processes in solids are discussed. The 
transition probabilities per unit time of the 3rd-order 
heavy particle assisted and the 4th-order electron as- 
sisted resonance-like double nuclear processes are deter- 
mined. The 3rd-order heavy particle assisted and the 
4th-order electron assisted resonance-like double nuclear 
processes may partly be responsible for the so called 
anomalous screening effect observed in low energy accel- 
erator physics investigating astrophysical factors of nu- 
clear reactions of low atomic numbers [l5|. The theo- 
retical description of the doubled processes discussed ex- 
tends the possible explanation and description of LENR 
with its nuclear physical background. It is found that the 
d + d — > 4 He process coupled to the f 6 Pd ^ x Ag 

process due to the quasi- resonant electron assisted dou- 
bled nuclear process has extremely large rate. The \He 
production with 34.7 MeV/He obtained in the leading, 
4th-order, quasi- resonant electron assisted d + d — > 2 ffe 
process fits well with the observed 32 ± 13 MeV/He value 
ofLENRs@. 

With the help of our theory describing resonance-like 
electron assisted doubled nuclear processes we estimated 
the nuclear power (|41|) created in the d + d — > 4 He pro- 
cess coupled to the ±§Pd + d — s- ^ 1 Ag process due to 
the quasi-resonant electron assisted doubled nuclear pro- 
cess and the nuclear power (|56|) due to double proton 
capture in Ni — H systems formed by Ni placed in H 2 
gas environment. The nuclear powers are consistent with 
observations. Moreover, the magnitude of the power ob- 
tained in Ni — H systems containing Ni in powdered 
form is considerable from practical point of view. 

The authors are indebted to K. Hartlein for his tech- 
nical assistance. 



where 



T 



(3) 



if 



EE 

\3 



yStyStyCb 



[E v — Ei) (Ep — Ei) 



(60) 



x l^L*(k 5 +k 6 ). 

ky denotes the wave vector of particle j in the final state. 
It will be seen that the process may have resonance like 
character if the masses of particles 5 and 6 differ signif- 
icantly, therefore the T^' attached to the leading graph 
(e.g. FIG. 1(a) in the case 7715 ^> m 6 discussed) is enough 
to calculate. kj n denotes the wave vector of particle j in 
state n — i, /i or v. The initial wave vector of particle 
2 is neglected in the Coulomb matrix element 

v cb _ 2 (2tt) 4 e 2 z 1 z 2 4>u (kig + k 2M ) 

V* 9s y3/2 |k 2 J 2 + A2 ' (bij 

where the 

&i(ki M + k 2M ) = 8^ 9 / 4 /3? /2 5(k lM + k 2M ) (62) 

formula can be used assuming that particle 1 is localized. 
tjjii stands for the Fourier transform of the initial state 
ipu- 

The two nuclear matrix elements 



V^ w (a) = -2/*/i4(M\/^ I 1 - - ) X (63) 



V 



V 



S (ki M - ki„) . 



and 



^(«) = - 2 / 2 / 23 (^)^ (l-l 



(64) 



are valid in the case of the process of FIG. 1(a) and 



(65) 



V^ W {b) = -2ff u {k lv )^j 



YItxR 
V 



2 

1- - 

e 



(66) 



APPENDIX I. RATE OF RESONANCE-LIKE 
HEAVY CHARGED PARTICLE ASSISTED 
DOUBLED NUCLEAR PROCESSES 



The rate of the process is 
2tt 



E 



(3) 



if 



S(E f - A) 



(59) 



stand for the process of FIG. 1(b). Here R is the ra- 
dius of a nucleon (we take R — 1.2 x 10 -13 [cm], that 
is the proton-radius) and the single nucleon approach in 
the Weisskopf approximation is used. The energy differ- 
ences in the denominator of (|60| are: the difference of 
the kinetic energies 



E,, — Ei 



2m\ 



2m 2 



Ei 



kin 



(67) 
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and 



{E v - E l ) a = 



ti 2 k 2v 

2m 2 



2 1,2 



2n%f 



A 5 



in the case of FIG. 1(a) and 



{E v - Ei) 



h 2 kl i h 2 kj 
2m\ 



2tok 



-E 



i.kin 



(68) 



(69) 



in the case of FIG. 1(b). In ([68]) and ([69]) the rest energy 
differences have to be also taken into account because 
of the nuclear reaction 1' + 4 — > 5 in the case of (a) and 
2' + 3 — » 6 in the case of (b) . A 5 = A i + A 4 — A 05 is the 
energy of reaction 1' + 4 — > 5, and A 6 = A 02 + A 3 — A 6 
is the energy of reaction 2' + 3 — > 6. The Aoj-s are 
the energy defects of the corresponding nuclei. The total 
reaction energy 



A s 



(70) 



After performing the E/ (really using the E/ 



dk.§ correspondence and 



J[V/(27r) 3 \ dk 5 xf[v/(2n) 3 _ 

integrating over k 6 ) in (|55|) the Dirac delta 5 (k 5 + kg) 
in (f6"0"]) will result k 6 = — k 5 , and the final energy can be 
written as Ef = h 2 k 2 / {2^0). The energy Dirac delta 
will result 



2^56A 

h 2 ' 



(71) 



Furthermore because of the presence of 8 (k 2 „ — k 2 „) 
and 8 (ki M — ki„) in the matrix elements (a) and 

V v f w (b) the £„,->/ [v/ {2nf] dk 2 „ and £„ -> 
/ F/(27r) d dki v will allow the k 

iv = k 2A , and ki 

v — ki^ 

substitutions in them and in the energy denominators 
as well. The initial kinetic energy Ei^in is neglected in 
{E v - E l ) a and {E u - E^. Thus 



{Eu - Ei) a = 



2m 2 



S a 



(72) 



with 



and 



8 a = ^A - A 5 
m 5 



{E v — E; 



h 2 k\ 

2m\ 



(74) 



one of 8 a and 8 is negative. Let us suppose that 8 a < 0. 
It means if 



fc2l,2 

211 = -8 a = A 5 - ^A 



2m 2 



m 5 



(77) 



then {E v — Ef) = and we find that the process (a) 



has resonance-like behavior when h 



2/i 



k r with k r 



2m 2 (^A 5 - -^AJ /h. If S a < then S b > and the 

process (b) can not have resonance character, and there- 
fore it is enough to calculate the Th ' attached to graph 
3(a). 

Let us introduce the half width T of the resonance with 
which the complex energy difference, denoted by suffix C, 



reads 



{E v - Ei) a)C 



2m 2 



S a -i- 



(78) 



that equals — if the resonance condition (1771) is 
met. With the use of the correspondence V, — > 

/ \V/ (2tt) 3 1 dk 2fl x J \v/ (2tt) 3 ] rfki M in (JBDJ and after 

carrying out the integral over ki^ the relation (f6"2"j) gives 
ki^ = — k 2A1 . Integrating over k 2AI we have the integral 
of the form 



I = 



h (flfyt) dk 2M 

9rr,„ ^ U Q 4 



(79) 



where /i (fc 2M ) is any function of fc 2/i , that in this case 

h (fc 2M ) = gs {k 2i , 2k 2ll ) / 23 (fc 2M )/i 4 (fc 2M ). (80) 

8 a is introduced 



For evaluating ([79]) the 5 (fc 2M 
and the substitution 



2m2 



1 



3 (fe/i) - 
is applied resulting 



P- 



1 



g{k2 ll ) 



iir8{g{k 2fi )) (81) 



(73) where 



P 



h (fc 2M ) rfk 2M 

fi 2 *- 2 

ft K2 M 1 X 
2m 2 T a 



(82) 



(83) 



and 



with 



Since 



8 b = **A 
m 6 



(75) 



(76) 



I s =tt j h {k 2fl ) 8 {g (fe/x)) rfk 2M . 
Using the identity 

<5(fc 2M -fc r ) 



(84) 



(85) 
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where k r (see ([9])) is the root of the g (fe/x) = equation 
(see ([77]) ). In a lower estimation of I as |7| = Ig, the 
relevant part of (|60|) is approximated as 



h (k 2 p) 



h 2 kl 
2m; 
-2;,r 



-dk. 



+ 5 a 



2ju 



(86) 



Applying the above relations and results, a lower approx- 
imation of the transition rate reads as 



Wf) = K^K[ 3 > (A, A 5 ) H (k r ) h 2 c 



(87) 



where K®>, K {3) (A, A 5 ) and H (k r ) are given by @, 
(JTJ) and ((SJ, and for h corr ^ and h corr ^ see dTUJ) and (TTT1) . 
that are defined by Eq.(45) of and fc r is determined 
by ©. 



VI. APPENDIX II. RATE OF 
RESONANCE-LIKE ELECTRON ASSISTED 
DOUBLED NUCLEAR PROCESSES 



The rate of the process 



with 



T^=Ti;'(a)+T^(b)+T^(c), (89) 



f£K 4) |V/-A) 



(88) 



-.(4)/ 



-.(4), 



where if/ 5 (a), T-^ (6) and ^(c) are the matrix ele- 
ments of the processes of FIG. 2 (a), (b) and (c), respec- 
tively. If processes of FIG. 2 (a) and (b) have resonance- 
like character (see later) then the process of FIG. 2 (c) 
has not and therefore its contribution may be neglected. 
Thus we take into account the contributions of 



T/St~\/Cb-\/St 

^ 4) («) = EEE 



pun 

(27T) 



(E p — Ei) (E v — Ei) 
<5(ki / + k 6 + k 7 ) 



(90) 



approximation (|62j) are used. Thus in each filled dot of 
the graphs representing a Coulomb interaction the mo- 
mentum (wave number) is conserved. The initial wave 
vectors of particles 1, 2 and 4 are neglected. The matrix 
elements of nuclear transitions 2', 3 — > 6 and 4', 5 — » 7 
are calculated in the Weisskopf approximation applying 
the appropriate one from formulae (j63[) . (1551) . (|6~4")l and 
(|66[) with the appropriate /23 (fc 2 ) and /45 (k 4 ) functions 
in it, respectively. The initial motion of particles 3 and 
5, i.e. their initial wave vectors are also neglected. In 
summing up for the intermediate and final states and for 
the square of the Dirac delta of argument of wave vector 
the correspondences and relations used above are applied 
again. (Remember, that now A 6 = A02+A03 — A 6 is the 
energy of reaction 2' + 3 — > 6 and A7 = A 5 + A 4 — A 7 
is the energy of reaction 4' + 5^7. The Acy-s are again 
the energy defects of the corresponding nuclei and the 
total reaction energy A = A$ + Aj). We calculate the 
rate of those processes in which the kinetic energy E\ / of 
the electron can be neglected in the energy Ef of the fi- 
nal state. Consequently the V/ (2ir) 3 J dkif will result a 
factor V/ (2-k) 3 (4tt/3) k\ Max , where k x . Max is the maxi- 
mum of the possible wave vectors of the electron in the 
final state. Furthermore, in the energy Dirac delta 



E f = E 6 + E 7 -A 



(92) 



is used. Neglecting also the final wave vector ki f of the 
electron in the Dirac delta, 5(k@ + k?) is used in (19"0"|) . 
(|9"Tj) resulting 



Ef 



h 2 k 2 7 

2fi 



(17 



(93) 



after integration over kg (ke = k 7 ). 

Since in the cases investigated ki/ is neglected, the 
wave number vector conservation in Coulomb scattering 
results the conservation of the magnitude of wave vector 
in lines 1', 2' and 4', and it is denoted by k±. Let us now 
investigate the energy denominators. The intermediate 
states are labeled with //, v and p. In the case of graph 
(a) 



(Ey — Ei 



and 



x/Stx/Stx/Cb 



p v n 

rCb 



{E p — Ei) (Eiy — Ei 



(Ep — Ei) V 



■<J (ki/ + ke + k 7 ) . 



(91) 



The outline of the calculation is the following. For the 
Coulomb matrix element of the process 1,2 — > 1', 2' we 
use the form V^ b given by Eq.(37) of [4]. For calcu- 
lating the matrix element of the Coulomb interaction of 
the process 1', 4 —> 1", 4' the form given by (|UTj) and the 



h 2 k 2 I 

(E v - E t ) a = ^-1 - A 6 + yf (hck 4 ) 2 + m 2 ^ - m e c 2 , 

(95) 



and 



{E P - E i ) a 



h 2 k 2 

2me 



h 2 k 2 
2m4 



(96) 



In obtaining E\ f is neglected. In the case of graph 
(b) 



(Ep - Ei) b - (Ep - Ei) a , 



(97) 
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(E v - Ei) b 



h 2 k 2 

2m 2 



H 2 kj 

2?7l4 ' 



and 



fc2p fc2p 



2ni4 



(98) 



(99) 



One can see from ((5BJ) and ((SSJ) that (£ p - £ l j) a = 
(£" p — Ei) b . Integrating over k7 and using the energy 
Dirac delta 



(E p - E t ,) a - (E p - Ei) b - 



h 2 k\ 
2m 4 



Sa,, 



with 



A 6 -^A 

m 6 



A fi - 



me + m 7 



(100) 



(101) 



If 5 a ,b > then in (j9TJ)) and (|9"Tj) resonance appears at fej 
given by (f2"5)) . Let us introduce again the half width T of 
the resonance with which the complex energy differences 
read 



h 2 k\ 



( E P - E *)a,C = ( E P ~ "iJb,C - 2mi 

The integration over k4 

f 23(h) J '45(h) 



S a . 



( io2 ) 



h 2 k 2 

2m 4 • 



-rfk4 



(103) 



= 4ir 2 k r 4 m 4 f 23 (kl)f u (k r 4 )/h 2 



with the aid of ([86|). where fc£ is given by (|25l) . 
(Now, because particle 1 is an electron gs{hv , h>i>) — 
gs(fa',h>) = 1 •) 

The energy differences in the denominators of (j9"0")) and 
(j9"T) will be 



(^-^) = <ya,6 (4/^^-1) =A , (104) 



(E v - Ek) b = —5 a ,b = A 6 , (105) 

A*24 



and 



<*-* ) ,_<^-^.-U£ + ,tesi)-A - . 



Now the rate of the 4th-order processes is 



(106) 



= TfWjb^ft^ftLr.e (107) 



where F G and x (A) are determined by <j23j), (j24|) 

and (|27p. fcj is given by (|23|) and for h corr ^ and h corr ^ 
see (Unj) and (JTTJ) . 
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